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ABSTRACT 

We present the first results from a very deep Chandra X-ray observation of the core of the 
Perseus cluster of galaxies. A pressure map reveals a clear thick band of high pressure around 
, the inner radio bubbles. The gas in the band must be expanding outward and the sharp front to 

^ ■ it is identified as a shock front, yet we see no temperature jump across it; indeed there is more 

\^ \ soft emission behind the shock than in front of it. We conclude that in this inner region either 

. thermal conduction operates efficiently or the co-existing relativistic plasma seen as the radio 

' mini-halo is mediating the shock. If common, isothermal shocks in cluster cores mean that we 

\ cannot diagnose the expansion speed of radio bubbles from temperature measurements alone. 

They can at times expand more rapidly than currently assumed without producing significant 
regions of hot gas. Bubbles may also be significantly more energetic. The pressure ripples 
found in earlier images are identified as isothermal sound waves. A simple estimate based on 
■ their amplitude confirms that they can be an effective distributed heat source able to balance 

Qh' radiative cooling. We see multiphase gas with about 10^ Mq at a temperature of about 0.5 keV. 

Much, but not all, of this X-ray emitting cooler gas is spatially associated with the optical 
^ • filamentary nebula around the central galaxy, NGC 1275. A residual cooling flow of about 

^ \ 50 M© yr^^ may be taking place. A channel is found in the pressure map along the path of 

• the bubbles, with indications found of outer bubbles. The channel connects in the S with a 

curious cold front. 
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1 INTRODUCTION 

The Perseus cluster, A 426, is the X-ray brightest cluster in the Sky 
and has therefore been well studied by all X-ray telescopes. The 
X-ray emission is due to thermal bremsstrahlung and line radiation 
from the hot IntraCluster Medium (ICM) and is sharply peaked on 
the cluster core, centred on the cD galaxy NGC 1275. Jets from 
the nucleus of that galaxy have inflated bubbles to the immediate N 
and S, displacing the ICM (Bohringer et al 1993; Fabian et al 2000). 
Ghost bubbles devoid of radio-emitting electrons, presumably from 
past activity, are seen to the NW and S. The radiative cooling time 
of the gas in the inner few tens of kpc is 2-3 hundred Myr, leading 
to a cooling flow of a few lOOMoyr^' if there is no balancing 
heat input. Energy from the bubbles or the bubble inflation process 
is a likely source of heat but the energy transport and dissipation 
mechanisms have been uncertain. 
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We have previously observed the Perseus cluster with the 
Chandra Observatory for 25 ks (Fabian et al 2000; Schmidt et al 
2002; Fabian et al 2002), 200 ks (Fabian et al 2003a,b, Sanders et 
al 2004, 2005) and now present here the first results from a further 
800 ks of observation. The total good exposure time is 900 ks. 

In the earlier work we discovered both cool gas and shocks 
surrounding the inner bubbles as well as quasi-circular ripples in 
the surrounding gas which we interpreted as sound waves generated 
by the cyclical bubbling of the central radio source. Related features 
have been seen in the Virgo cluster (Forman et al 2003). The NW 
ghost bubble has a horseshoe-shaped optical Ha filament trailing it 
which we interpret as showing the streamlines in the ICM. On this 
basis we concluded that the ICM is not highly turbulent and thus 
that viscosity is high enough to dissipate the energy carried by the 
sound waves (Fabian et al 2003a,b). Such an energy transport and 
dissipation mechanism is roughly isotropic and can thereby provide 
the required gently distributed heat source required by observation 
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Figure 1. Colour image assembled from separate images in the 0.3-1.2 
(red), 1.2-2 (green) and 2-7 keV (blue) bands. 



of this and other similarly X-ray peaked clusters (Ruszkowski et al 
2004a,b,2005; Reynolds et al 2005; Fabian et al 2005). 

Our goal in the present work is to determine the temperature 
and pressure of the ICM accurately so that we can study the pro- 
cesses taking place there in more detail. We indeed confirm that the 
pressure jumps at the weak shock surrounding the inner bubbles 
and also that the ripples represent significant ripples in pressure. 
The temperature does not jump at the shock, however, which may 
be due to the action of efficient thermal conduction. The energy 
from the bubbles propagates through isothermal sound waves and 
conduction in the inner regions. If this is a common property of 
such regions then some of the otherwise puzzling behaviour can be 
understood. 

The redshift of the Perseus cluster is 0.0183, which for a Hub- 
ble constant of 71 km s^' Mpc^' corresponds to a luminosity dis- 
tance of 78.4 Mpc and an angular scale of 367 pc per arcsec. 



2 THE DATA 

The Chandra observations used for the analysis presented in this 
paper are listed in Table Q The total exposure time of just over 
1 Ms is reduced after removing periods containing flares to 890 ks. 
To filter the datasets we examined the lightcurve between 2.5 and 
7 keV on the ACIS-Sl CCD. The SI CCD is back-illuminated like 
the S3 and so it is best CCD to search for flares, as the Perseus 
cluster emission is dominant over flares on the S3 CCD. The CIAO 
3.3.2 LC_CLEAN tool was used to remove periods from the dataset 
which deviated away from the median count rate of all the obser- 
vations. Observations 3209 and 4289 did not use the S 1 CCD and 
did not show any flares on the S3 CCD, and so were left unfiltered. 

The level 1 event files were reprocessed using the PSU CTI 
corrector (Charge Transfer Inefficiency; Townsley et al 2002a, 
2002b). Level 2 event files were produced by removing standard 
grades and bad time intervals. Each of the event files was then re- 
projected to match the coordinates of the 04952 observation. Im- 
ages of the data in this paper were produced by summing all the 
images from the individual datasets. To correct for exposure vari- 
ation we created exposure maps for each of the CCDs for each of 



the datasets and for each of the bands. The summed images were 
then divided by the summed exposure maps. 

We have produced unsharp-mask images by subtracting im- 
ages which have been smoothed on two lengthscales. Fig. 2 (top) 
shows the result after using Gaussian smoothing of 2 and 20 pix- 
els. The ripples are very clear, out to radii of 3^ arcmin (60- 
80 kpc) from the nucleus. An arclike step in surface brightness oc- 
curs ~ 1.5 arcmin S of the nucleus. A cold front is seen to the SE 
(we verify that the pressure is approximately continuous across the 
sharp surface brightness change in Section 4). Such features were 
first seen in Chandra images of clusters by Markevich et al (2000). 
There is a major difference with the feature seen here, however, 
since it is concave and cannot be due to the core moving through 
a wider hotter gas. It does however appear to connect the 'bay' to 
the S of the nucleus which connects in towards the nucleus along a 
narrow channel emerging to the SSW from the inner regions. This 
corresponds spatially to a weak outer Ha filament (Sanders et al 
2005) although extends much further than any optical emission. X- 
ray emission is also associated with a much more dominant long 
radial Ha filament seen to the N of the nucleus (see e.g. Conselice 
et al 2001). The X-ray feature appears to break beyond the ripples 
and is labelled as the Ha fountain. 

As will be discussed later, we suspect that the radial features 
are due to cold and cooler gas dragged out from the centre by rising 
buoyant bubbles. They represent the main axis along which most 
of the bubbles rise. The S cold front could then be the edge of a 
giant hotter bubble either produced by a past major outburst of the 
nucleus (cf. McNamara et al 2005 for Hydra A) or where the hot 
gas accumulates due to the interior entropy of the bubbles matching 
the external value there. 



3 TEMPERATURE AND PRESSURE MAPS 

The total of ~ 70 million counts in the final all band image from 
the ACIS-S3 chip means that we can measure spectral properties 
on unprecedented small scales. In order to proceed we have divided 
the image into bins with approximately the same number of counts 
and used XSPEC 11.3.2 (Arnaud 1996) fitting with MEKAL models 
(Mewe, Gronenschild & van den Oord 1985; Liedahl, Osterheld & 
Goldstein 1995) to obtain spectral parameters, fitting between 0.5 
and 7 keV. The temperature map shown in Fig. |4] was derived in 
this way using a contour binning approach (Sanders in preparation) 
with 625 or greater ct per spectrum. In each fit the metallicity (in 
Solar ratios; Anders & Grevesse 1989) and absorption column den- 
sity were fixed at values measured when fitting spectra from bins 
containing 10'* ct per bin or greater, except in the region around the 
High Velocity System where the absorbing column density was al- 
lowed to be free. The results for these parameters is broadly similar 
to those found in our earlier work (Sanders et al 2004). Details will 
be given in a later paper. For the present work we concentrate on 
the temperature and emission measure distributions. 

We used standard blank sky observations to act as back- 
grounds for the spectral fitting. The background observations were 
split into sections to match the ratio of exposure time between 
each foreground observation. These datasets were then reprojected 
to match the foreground observations, and then reprojected to the 
4952 observation. The exposure time of the backgrounds were al- 
tered to ensure the same rate of counts between 9 and 12 keV as 
their respective foregrounds, in order to correct for the variation 
of background with time. To create a total spectrum, the spectra 
from each of the individual observations were added together, ex- 
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Obs. ID 


Sequence 


Observation date 


Exposure (ks) 


Nominal roll (deg) 


Pointing RA 


Pointing Dec 


3209 


800209 


2002-08-08 


95.8 


101.2 


3:19:46.86 


-1-41:31:51.3 


4289 


800209 


2002-08-10 


95.4 


101.2 


3:19:46.86 


-1-41:31:51.3 


6139 


800397 


2004-10-04 


51.6 


125.9 


3:19:45.54 


-1-41:31:33.9 


4946 


800397 


2004-10-06 


22.7 


127.2 


3:19:45.44 


-1-41:31:33.2 


4948 


800398 


2004-10-09 


107.5 


128.9 


3:19:44.75 


-1-41:31:40.1 


4947 


800397 


2004-10-11 


28.7 


130.6 


3:19:45.17 


-1-41:31:31.3 


4949 


800398 


2004-10-12 


28.8 


130.9 


3:19:44.57 


-1-41:31:38.7 


4950 


800399 


2004-10-12 


73.4 


131.1 


3:19:43.97 


-1-41:31:46.1 


4952 


800400 


2004-10-14 


143.2 


132.6 


3:19:43.22 


-1-41:31:52.2 


4951 


800399 


2004-10-17 


91.4 


135.2 


3:19:43.57 


-1-41:31:42.6 


4953 


800400 


2004-10-18 


29.3 


136.2 


3:19:42.83 


-1-41:31:48.5 


6145 


800397 


2004-10-19 


83.1 


137.7 


3:19:44.66 


-1-41:31:26.7 


6146 


800398 


2004-10-20 


39.2 


138.7 


3:19:43.92 


-1-41:31:32.7 



Table 1. The Chandra observations included in this The exposure given is the time remaining after filtering the lightcurve for flares. The observations were 
taken with the aimpoint on the ACIS-S3 CCD. Positions are in J2000 coordinates. 




Figure 2. Unsharp mask image made from the whole 0.3-7 keV band by subtracting an image smoothed with a Gaussian of dispersion 10 arcsec from one 
smoothed by 1 arcsec and dividing by the sum of the two images. Various features are labelled on the lower contrast image at the left. 



eluding observations which did not have any counts in the region 
examined. The background spectra were added together similarly. 
The standard PSU CTI corrector response was used. Ancillary re- 
sponses for each dataset and region were produced using the CIAO 
MKWARF tool, weighting CCD regions using the number of counts 
between 0.5 and 7 keV. These ancillary responses were averaged 
for each region, weighting according to the number of counts be- 
tween 0.5 and 7 keV for a particular dataset. 

The temperature map (Fig. 4) shows in great detail the 'swirl' 
around NGC 1275 (Churazov et al 2000). Whether the swirl is re- 
ally a single connected structure or an outer ring partially opened 
on the E and connected to the rim of the inner N bubble (Dunn et al 
2005b) remains unclear. Some 'fountaining' can be seen to the N 
of this N bubble. This is associated with the N optical Ha filaments 
which are surrounded by gas at about 1 keV (see Fig. 3). A disrup- 
tion in the outer ring is seen to the SE of the nucleus coincident 
with the optical 'blue loop' ; this is discussed in Section 6. 

We now focus on measurements of entropy S and particularly 
the pressure P of the gas. A simple method for obtaining these 



quantities is to assume that the density n is proportional to the 
square root of the X-ray surface brightness and then use P = nkT 
and S = Trc'^l^ . Here we use a slightly better approach based on 
the emission measure. A, obtained from spectral fits. This is pro- 
portional to rP'V where V is the volume along the line of sight. 
Since the emission is strongly peaked we ignore V at this stage and 
produced 'projected' entropy and pressure maps (Fig. 5). 

The entropy map (Fig. |5| left) emphasises where gas may 
have cooled and resembles the temperature map. The pressure map 
(Fig.|5]right) on the other hand shows clearly a thick band around 
the inner radio bubbles and little sign of azimuthal asymmetry. As 
found by Sanders et al (2004) the pressure distribution is reason- 
ably circularly symmetric, as expected for gas close to hydrostatic 
equilibrium. This is not just a consequence of our volume assump- 
tion, since we see that the 'swirl' in temperature has completely 
disappeared, as has the arc noted in Fig. 2. 

A thick, higher pressure band surrounds the radio-filled cavi- 
ties or bubbles (Fig. 6). This presumably is shocked gas produced 
by the inflation of the bubbles. It is remarkable that we see it as two 
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Figure 3. Colour image made from the 0.3-1.2 (red), 1.2-2. (green), 2.-7 keV (blue) bands. A 10 arcsec smoothed image has been scaled to 80 per cent of its 
intensity and then subtracted in order to bring out fainter features lost in the high intensity range of raw images. The blue structure to the N of the nucleus is 
caused by absorption in the infalling high velocity system, projected at least 60 kpc in front of the nucleus of NGC 1275 (Gillmon et al 2004). 



mostly complete rings in the projected pressure map. This means 
that the two bubbles cannot lie in the plane of the Sky but must be 
arranged so that one is nearer us than the other. Since the nearer 
radio jet is the S one based on VLBI radio data, we suppose that 
the nearer bubble is the S one. 

There is some azimuthal asymmetry in the pressure map, 
mostly associated with the bubbles. In order to see this we have 
subtracted the mean pressure at each radius to produce the pressure 
difference map (Figs |6| and 0. There are clearly some lower pres- 
sure regions to the N and S, probably associated with older, outer 
bubbles. The region to the SSW has a higher metallicity likely due 
to older bubbles dragging metal-rich gas there (Sanders et al 2005). 

To the south we see two further tangential arclike pressure 
minima beyond the outer S bubble. These coincide with the high 
abundance shell reported by Sanders et al (2005). To the N we also 
see a large arclike pressure minimum. 

We suspect that these arclike pressure minima are old bubbles. 



The large size of these bubbles could indicate that the activity was 
much stronger in the past, so blowing larger bubbles, or may just be 
due to bubbles merging. Surrounding gas may leak into the bubbles 
so making them less buoyant, or magnetic fields and structures may 
be important. 

No pressure jump is associated with the concave structure to 
the S, confirming that it is part of a cold front. 



4 THE SHOCK AND RIPPLES 

We now discuss the detailed behaviour of the temperature, density 
and pressure around the shock surrounding the inner bubbles and 
the ripples. Sectors have been defined to the NE, E, S and NW of 
the nucleus and spectra extracted from bins spaced 5.4 arcsec in 
radius (Table|2j- 

The projected temperature profiles are shown in the top pan- 
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Projected Entropy Projected Pressure 

Figure 5. Entropy (left) and pressure (right) maps. The entropy map was calculated using kTA-^/^, in units of keV cm^/^ arcsec^''', where A is the MEKAL 
normalisation per square arcsecond. The pressure map was calculated using kT A^^^, in units of keV cm^'^^ arcsec^'. These maps were generated by fitting 
regions containing approximately 625 counts or greater. 



Name Centre RA Centre Dec Start angle Stop angle 



Ripples in pressure ranging from ±5 to ±10 per cent in pressure 
are seen out to 50 kpc or more. Such pressure variations cannot 
be static. They resemble sound waves so, following our earlier 
work (Fabian et al 2003a), we interpret them as sound waves' pro- 
duced by the cyclic bubbling behaviour, or at least major expansion 
episodes, of the inner bubbles. In the inner region they are high 
pressure regions fronted by a weak shock, further out the shocks 
weaken and are not distinguished from the overall pressure distur- 
bance or ripple. 

A simple calculation serves to show the potential for the rip- 
ples to heat the gas provided that viscosity is high enough to dissi- 
pate their energy. Let us consider the region within 50 kpc where 
the ripples are most clearly seen. If the ripples move at 1000 km s^ ' 
then they cross this region in 5 x 10^ yr. They cause the pressure to 
oscillate with an amplitude of 5-10 per cent, which we conserva- 
tively take as 5 per cent of the thermal energy there. Consequently 
they can balance cooling provided that the cooling time (roughly 
for the gas to lose all its thermal energy) is 20 times the crossing 
time or a Gyr. This condition is well met since the cooling time of 
the gas at 50 kpc is 2-3 Gyr dropping to about 2 x 10^ yr near the 
centre (for the hotter gas). The waves therefore need to dissipate 
about half their energy by 50 kpc. 

The vertical dashed line in the upper panel of Fig. 8 corre- 
sponds to the abrupt edge to the thick pressure band around the N 
bubble. It is most clearly seen in this direction because there is little 
lower energy emission superimposed upon it as appears to happen 

' Or at least they operate as sound waves in a mixed thermal/relativistic 
plasma 



North-east 03:19:48.11 +41:30:41.22 22.5 52.9 

East 03:19:48.11 +41:30:41.22 91 106 

South 03:19:45.92 +41:30:19.58 136.9 164.2 

North-west 03:19:48.11 +41:30:41.22 294.4 334.2 

Table 2. Sectors used to generate temperature, density and pressure pro- 
files. Angles are measured from North in the eastern direction. Coordinates 
shown are J2000. 



els in Fig. 8. We see generally that, apart from where the inner and 
ghost bubbles lie, the temperature profiles are smoothly increasing 
from the inside out. We have explored deprojected temperatures us- 
ing the XSPEC PROJCT routine but find them unstable with the tem- 
peratures depending on which bins are used. Dropping bins from 
either the outside or the inside can have large unpredictable effects 
on the results found for bins at intermediate radii. This may be due 
to the real geometry being obviously different from the spherical 
or ellipsoidal geometry assumed by the routine or to the gas being 
multiphase. We therefore proceed to combine the projected temper- 
atures with approximately deprojected densities from the emission 
measures to obtain pressure profiles. The deprojected densities are 
calculated by subtracting the contribution to the fitted normalisa- 
tion at each radius by that expected from outer shells, assuming 
projection in a spherical geometry. 

The density and thus pressure profiles show variations (Fig. 8) 
corresponding to the ripples seen in the unsharp mask image 
(Fig. 2). The pressure residuals from a smooth power-law fitted 
from 20 to 70 kpc are shown in the lowest panel of the plots, to- 
gether with the residuals predicted from the unsharp mask image. 
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Figure 8. Temperature, density, pressure and pressure variation profiles. Tlie red line shows the ripples from the unsharp mask image. The top left figure shows 
the profile in the NE direction, top right shows E direction, bottom left shows S direction and bottom right shows NW direction. The dashed line in the NE 
profiles indicates the position of the shock front. 
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Figure 4. Temperature map calculated by fitting spectra with approximately 
625 counts or greater. The uncertainties of the individual fits range from 
8 per cent in the coolest regions to 20 per cent in the hottest parts, ignoring 
the uncertainty on the metallicity and absorbing column density. 
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Figure 7. Thermal pressure map where the mean pressure at each radius has 
been subtracted. In this map the temperatures and normalisations were mea- 
sured using regions containing approximately 10* counts or greater Note 
that the 'channel' caused by a sequence of 4 thermal pressure dips running 
to the S of the nucleus. The outer ones are assumed to old ghost bubbles and 
the missing pressure is assumed to be due to relativistic plasma. A twisted 
channel is also seen to the N. 





Figure 6. 1.4 GHz Radio map in blue superimposed on the pressure dif- 
ference map in red, where the average pressure at each radius has been 
subtracted. In this map the temperatures and normalisations were measured 
using regions containing approximately 10* counts or greater. 



at other azimuths. Our results show it to be associated with a jump 
in density and this pressure and should be a shock front. However 
the projected temperature hardly changes, even dropping shghtly 
postshock (Fig. 8), whereas it should rise on the basis of the factor 
1.39 density jump, from 4.04 keV to about 5.1 keV if the gas is 
adiabatic (see also Fabian et al 2003a). 

The very deep image enables us to now examine the exact tem- 
perature structure across this shock front. Because of the geometri- 
cal restrictions to any extended deprojection approach, mentioned 
above, we compare spectra either side of the shock. All methods, 
including single and two temperature, or using the outer emission 
as background for the inner, show no sign of any hotter temperature 
component within the shock in the range of 4—6 keV. The ratio of 
the spectra post and preshock, plotted in Fig. 9, clearly shows that 
the emission inside the shock is softer and inconsistent with the 
predicted rise, if y = 5/3. All the spectral fits have been carried out 
assuming that the gas is in coUisional and ionization equilibrium 
and homogeneous within each region. 

We have also carried out a multi-temperature analysis of the 
emission both inside and outside the shock (Fig. 10). 56 per cent of 
the counts per unit area across the jump going from out to inside 
is in an increased 4 keV component and 28 per cent in a 2 keV 
one. There is good evidence that the gas is multiphase. The 4 keV 
component, which is the volume-filling one, shows an increase of 
density across the shock but no increase in temperature. 

The soft X-ray emission from the outer radial, N optical Ha 
filament stops abruptly at the shock front. It is possible that the 
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Figure 9. Top: The ratio of the spectrum inside the shock to that outside. 

The spectra have been normalised to have the same number of counts. The 
sohd hne shows the ratio is the abrupt density jump corresponds to a weak 
shock with 7=5/3. Bottom: Regions used in the above spectral analysis, 
shown superimposed on an unsharp-masked image. The lower region is that 
within the shock, the upper one is outside it. 



2 keV component arises from shock heating of this cooler (1 keV) 
gas. The pressure jump at the shock front indicates that this is a 
weak shock so most of the heating is just PAV compression. 

Why the 4 keV component does not change is puzzling if the 
gas is adiabatic; much of the heating in a weak shock is just com- 
pressional PAV heating. 

A fair approximation for the temperature behaviour expected 
in a shock where the density jumps by the observed factor of 
1.39 (Fig. 8) from 4.04 keV is that the postshock temperature 
r = 2.5 + 1.567keV. The observed behaviour of the gas is there- 
fore explained if the 4 kcV component is isothermal with 7~ 1. 

A simple explanation for this would be that thermal conduc- 
tion is operating on this volume-filling gas phase. The electrons. 







o 


Inner annulus 


0.01 




A 


Outer annulus 



6 
o 
Z 



10- 



10' 



10- 



CD 



® 



1 

Temperature (keV) 



10 



Figure 10. Distribution of temperature components within and outside the 
shock, using the same regions as shown in Fig. 9. The regions have been 
fitted with models having temperature components at 0.25, 0.5, 1. 2. 4. 8 
and 16 keV. Note that both the 4 and 2 keV components increase within 

the shock, with the 2 keV one increasing proportionately more. Only upper 
limits are obtained for the 0.25, 1 and 8 keV components (the ends of the 
error bars are shown as the vertical Unes). 



moving faster than the ions, can go ahead of the shock (see e.g. 
Zel'dovich & Raizer 1966; Borkowski, ShuU & McKee 1989). If 
the magnetic field in this region is mostly radial, then conduction 
can eliminate temperature differences on a timescale 



fcond ~ nkf/K = 2.3 X lO^n^^r-^/^ yr 



(1) 



which compares with the timescale for distance t of matter to ac- 
cumulate behind the shock (in the rest frame of that matter) 



<shock = 4.8xl0«'£yr 



(2) 



where I is the length inside the shock in units of 3 kpc (the bins 
in Figs. 8 and 9 are 2 kpc apart). The full Spitzer (1956) rate for 
conduction is assumed here and the total postshock density is in 
units of the observed value of 1 cm^^ (Fig. 8), the postshock tem- 
perature (assuming an adiabatic gas) is in units of the expected one 
of r = 5.1 keV = 6.3 x 10^ K. Magnetically-isolated blobs such as 
may comprise the lower temperature component are compressed 
adiabatically. The post-shock timescale for electron-ion equilibra- 
tion is comparable to the above time (~ 2 x 10* yr). This enhances 
the effect of conduction (halving it) since only the electron energy 
needs to be conducted ahead of the shock. We envisage that the ion 
temperature jumps at the shock front but that the electron tempera- 
ture varies smoothly through this region, with a hotter precursor ex- 
tending into the unshocked gas (see Fig. 7. 19 on p5 19 of Zel'dovich 
& Raizer 1966); both the electron and ion densities jump at the 
shock front. 

This result introduces the possibility that thermal conduction 
is effective in parts of the innermost regions of clusters. It has been 
proposed and tested as a means for heating the gas from the outside, 
but found to be inadequate for clusters and regions below 5 keV 
(Kim & Narayan 2003; Voigt & Fabian 2004). What is needed in 
the Perseus cluster is for thermal conduction to operate throughout 
much of the inner hotter volume-filling phase. The ripples would 
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therefore be isothermal sound waves (see Fabian et al 2005 for a 
comment on this possibility). Both sound waves and conduction 
are then effective in distributing the PdV work done by the bubbles 
into the surrounding gas. Repeated bubbling in the central region 
may have ordered the magnetic field into a roughly radial structure. 

Cooler and/or cooling temperature components embedded in 
the hotter gas behind the shock can damp any temperature rise be- 
hind the shock if they mix with the hotter gas. The mass fraction 
of cooler gas required (approx 30 per cent) appears not to be high 
enough (see Fig. 10) for this process to be important. It remains 
possible that mixing takes place with larger masses of unseen cold 
gas which radiates much of the thermal energy in yet imseen bands. 

An issue which could be very important for shock propaga- 
tion in the irmer intracluster medium is the presence of a relativis- 
tic plasma (cosmic rays and magnetic field) in the inner core of the 
Perseus cluster. This is evident here from the synchrotron emission 
seen as the radio 'mini-halo' (Pedlar et al 1990; Gitti et al 2003) 
and the inverse Compton emission seen as a hard X-ray flux com- 
ponent (Sanders et al 2005; it appears as the 16 keV component in 
Fig. 10). In the collisionless conditions relevant to the shock it may 
be possible that the relativistic plasma soaks up the energy, leav- 
ing the gas isothermal. Indeed it could be repeated shocks from the 
bubbles which reaccelerates the relativistic particles. They could re- 
distribute the energy to larger radii, serving to transport some of the 
energy and creating a distributed heat source for the gas. We note 
that the electron temperature observed behind the strong shocks in 
young supernova remnants do not always fit expectations for sim- 
ple hydrodynamical shocks, probably due to particle acceleration 
(Rakowski 2005 and references therein). Although promising as a 
mechanism, there are many uncertainties as to how it could operate 
and why there is no sharp rise in either the synchrotron emission 
seen in radio maps nor in the inverse Compton emission at the po- 
sition of the shock. Moreover it does not explain how the electrons 
avoid compressional heating. 

The isothermal nature of the inner gas raises the possibil- 
ity that the bubbles expand much faster than previously suspected 
from observations. Initial models for the action of a central radio 
source on the ICM by Heinz et al (1998) predicted that the bubbles 
would be surrounded by shocks but Chandra showed no evidence 
for shock-heated gas. Efficient thermal conduction will however 
eliminate shock heating as a diagnostic. Consequently the bubbles 
may expand, at times, faster than inferred, even supersonically. The 
likely behaviour given the variability of radio sources is that they 
expand in fits and starts, with each rapid expansion phase giving 
rise to a sound wave. The observation from the Perseus cluster of 
only one set of ghost bubbles within 50 kpc radius yet 3 or more 
ripples allows for each bubbles to generate several ripples before 
becoming buoyant enough to separate and rise. This means that 
any estimate of bubbling power based simply on buoyancy times 
(e.g. Birzan et al 2004; Dunn et al 2004, 2005a) is a lower limit. 

A further issue with regard to the energy injected by the bub- 
bles is the thickness of the postshock gas. This is very similar to 
the radius of the bubbles, so has a volume about 7 times that of the 
bubbles themselves. The pressure in the shocked gas is 30 per cent 
above the outer unshocked gas, P, (Fig. 8) so the energy content of 
the postshock gas is more than twice that obtained by assuming it is 
just PV where V is the volume of a bubble. The work done (/ PdV) 
will be yet higher if some has been transported away by conduction 
or relativistic particles. 




1 2 
Tempemture (keV) 

Figure 11. Distribution of mass at the fixed temperatures from within the 
innermost 1.5 arcmin radius. For comparison the expected result from a 
constant pressure cooling flow is shown. 



5 THE MULTIPHASE NATURE OF THE GAS 

Figs. 1 and 3 clearly show filamentary soft X-ray emission which 
is closely associated with the optical Ha filaments (Fabian et al 
2003b). This soft emission has a temperature of between 0.5 and 
1 keV and would appear much brighter if the Galactic column den- 
sity to the Perseus cluster were not as high as the observed value of 
~ 1.3 X 10^' cm~-^ . The Doppler velocities determined for the fila- 
ments are 100 — 200 km s~^ and coherent over many kpc (Hatch 
et al 2005b) so, given their large radial extent and likely origin 
as being pulled out from the centre by rising bubbles (Fabian et 
al 2003b), the lifetimes of the filaments are several tens million 
yr, or even longer. In order to survive in the surroimding hot gas 
they must be insulated from it or thermal evaporation would have 
caused them to disappear within a million yr (equation 1). Con- 
duction must therefore be highly suppressed, by at least a factor of 
100, probably due to magnetic fields along their length (conduction 
is suppressed perpendicular to the field direction). 

As already mentioned, the filaments stop at the shock which 
is probably disrupting them there. The filaments coincident with 
the shock to the SE are probably just projected in front of the 
shock and are not within it. Such magnetically isolated regions need 
not completely vanish once they are disrupted and may survive as 
higher temperature blobs maintaining their isolation. The gas can 
therefore be multiphase, not due to a thermal cooling instability, 
but to the forced mixing of different components. Whether there 
is then slow conductive evaporation or radiative condensation (see 
e.g. Bohringer & Fabian 1989) or turbulent mixing (e.g. Begelman 

6 Fabian 1990; Loewenstein & Fabian 1990) remains to be seen. 

We have therefore conducted a multi-temperature determina- 
tion of the gas distribution in the Perseus core. The individual spec- 
tra generated from regions chosen to contain 10'* counts or greater 
have been fitted with a multi-tcmpcraturc model consisting of gas at 
0.5, 1, 2, 4, 8 and 16 keV (see also Sanders et al 2004 for similar fits 
to the 200 ks data). The results have been mapped in terms of mass. 
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Figure 12. The Ha filaments from the optical narrow-band image of Conselice et al (2000) is shown for comparison with maps of mass in the 0.5, 1, 2 and 
4 keV components. 



determined from the emission measure n^V divided by the density 
n relevant for the pressure at that radius and measure temperature 
of the component. Of course the volume filling factor of the gas 
significantly different from the mean temperature found from the 
single temperature fits (Fig. 4) is small. The widely differing mass 
distributions on the Sky show that the gas is genuinely multiphase 
(i.e. having different temperatures at the same radius) and we are 
not mapping a mere projection effect. 

We see a striking similarity in the 0.5 keV map to the op- 
tical filaments with a total mass in the hot gas much larger than 
typically found from an estimate of 3 x 10^ Mq, based on the total 
luminosity of Ha (Heckman et al 1989), a temperature of 5000 K 
for the gas smd the surrounding pressure found here for the outer 
filaments. The continuing pressure rise to the centre will reduce 
this estimate and addition of molecular hydrogen, seen in the in- 
frared (e.g. Hatch et al 2005b and references therein) will tend to 
increase it, so it should be a reasonable estimate. The mass maps at 



the various temperatures are plotted in Fig. 12 and the total mass 
distribution, determined from the inasses within 1.5 arcmin of the 
nucleus, is shown in Fig. 11. The different points at each tempera- 
ture show the total mass including only those which are significant 
to 1, 2 and 3-<t. Noise will be a strong contaminant to the lowest 
significance point. For comparison, the mass distribution expected 
for a steady cooling flow of 300 Mq yr^^ is superimposed. 

Interestingly, we find that there is a large drop-off in mass at 
1 keV but a recovery at around 0.5 keV. This rise is of course due 
to the filament region. Until we know the fate of such material, 
in terms of whether it is being heated or cooled by radiation or 
mixing, we cannot say whether the bulk of the cooler gas, which 
lies in an E-W extended clump around the nucleus, is the residual 
of a cooling flow or not. We note that Bregman et al (2005) find OVI 
emission (characteristic of gas at 5 x 10^ K) in a 30 arcsec FUSE 
aperture centred on the nucleus consistent with a mass cooling rate 
of about 50 Mq yr~^ . This is comparable to the rate inferred from 
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Figure 13. The blue loop is shown in X-ray temperature (left), Ha (from 
ConseUce et al 2001, centre) and blue Ught (from JKT; right). 

our mass determination from gas at 0.5 keV (i.e. ~ 5 x 10^ K), since 
the radiative cooling time of gas between 0.5 and 2 keV in the iimer 
parts of the cluster is about 10** yr, where T is in keV. Peterson 
(2005; private communication) finds a limit of only 2OM0yr^' 
from a search for FeXVII emission in XMM-Newton RGS spectra 
of the iimer 30 arcsec radius. 

The fact that we see less gas at 1 keV could be the conse- 
quence of cooling due to mixing, rather than radiation, dominating 
in that temperature range. Such a possibility has been discussed by 
Fabian et al (2002) and Soker, Blanton & Sarazin (2002). The en- 
ergy of the hotter gas could in part go to heating the cooler gas at 
~ 10** K where there has long been a heating and excitation prob- 
lem (Johnstone, Fabian & Nulsen 1987; Heckman et al 1989; Sabra 
et al 2002). Indeed a mixing solution similar to a turbulent, radia- 
tive mixing layer seems inevitable given the much lower mass in 
cold gas below 10"* K then at 0.5 keV. 

A final inter-relationship between the hotter X-ray emitting 
gas and the optical filaments is shown in Fig. 13. There is a par- 
tial ring structure to the SE in the temperature swirl, resembling a 
letter 'C written backwards. It coincides with some bright optical 
filaments and in particularly with the 'blue loop', first remarked on 
by Sandage (1972) and seen well in many recent images (e.g. the 
blue band Jacobus Kapteyn Telescope, JKT, image of Fig. 12). We 
presume that gas in the swirl at this location collapsed and formed 
the stars in the astronomically recent past. 

Heinz & Churazov (2005) have proposed that the relativistic 
component discussed in Section 4 could exist in small blobs which 
could help to dissipate wound waves. We see no obvious signs of 
small holes in the X-ray emission larger than a few 100 pc in size. 
How well the relativistic and thermal components are mixed is of 
importance for transport processes in the region. 



6 DISCUSSION 

We have found that the shock seen in our 200 ks image is isother- 
mal. The ripples seen beyond the shock are therefore likely to be 
isothermal waves. Their energy is then dissipated by viscosity. Con- 
duction and sound waves can act together to dissipate and distribute 
the energy from the radio source, and ultimately the central massive 
accreting black hole. An isothermal shock allows energetic bub- 
bling to occur at the centre without overheating the iimermost re- 
gion, a problem noted by Fujita & Suzaki (2005) and Mathews et 
al (2005). 

In the work of Fujita & Suzaki (2005) it is assumed that all 
wave dissipation occurs at the shock front and does not include any 



later dissipation via viscosity as the (observed) waves propagate 
further. In one model they include conduction at 20 per cent of the 
Spitzer rate and find agreement with the shape of the temperature 
and density profile. However most of the energy in their model is 
supplied by thermal conduction from the hotter outer gas; the AGN 
only dominates over the region from 20-30 kpc where the shock 
occur. As they remark, a double heating model with the AGN heat- 
ing the inner regions and conduction the outer was proposed earlier 
by Ruszkowski & Begelman (2002). 

It is clear from the temperature profile shown in the top left 
panel of Fig. 8 that conduction of heat from the outer hotter gas is 
not important within the inner 60 kpc in the Perseus cluster, since 
the temperature profile is so flat. Indeed from 40-55 kpc the gradi- 
ent acts in exactly the wrong direction. As discussed in Section 
4, the observed ripples (which are strong sound waves or weak 
shocks) have more than sufficient energy to heat the inner 50 kpc 
and so it is not clear that any thermal conduction of heat from the 
outer gas is required. What our analysis has shown is that thermal 
conduction acting in the inner regions can account for the observed 
isothermal nature of the shock and so prevent the problem of an ac- 
cumulation of hot shocked gas. The conduction merely acts to me- 
diate the shock and redistribute the energy from the central AGN. 

The magnetic configuration of the field in the core is crucial to 
the conductive behaviour. We require an approximately radial field 
across the shock, which is not understood. One possibility is that 
it arises as a consequence of cooling and compression of the inner 
gas in the past which leads to the frozen-in magnetic field being 
predominately radial (Soker & Sarazin 1990). Nearby we have low 
temperature Ha emitting filaments which must be many 10s mil- 
lion yr old and so magnetically isolated. We also find evidence for 
multi-temperature, presumably multiphase, gas. The magnetic con- 
nectedness is crucial to how the gas behaves. It raises the possibility 
that the swirl seen in the temperature and entropy maps is magneti- 
cally separate from the rest of the gas. Perhaps it is a fossil from the 
merger of a galaxy with NGC 1275 where the incoming (less mas- 
sive) galaxy 'combed' the field into the apparent swirl. The 2 keV 
gas immediately around the rim of the N inner bubble is presum- 
ably protected from evaporation by a tangential field there. Such 
speculation may eventually be testable when it is possible with the 
Expanded Very Large Array to carry out Faraday Rotation studies 
in this region at higher frequencies and greater sensitivity than cur- 
rently feasible. Preliminary indications from high resolution stud- 
ies of the nucleus with the Very Long Baseline Array indicate fairly 
extreme Rotation Measures of up to 7000 radian w? (Taylor et al. 
2005, in preparation). 

We have also found a roughly N-S channel in the pressure dif- 
ference map which demonstrates the passage of a sequence of radio 
bubbles. The outer ones are large and could be where they accumu- 
late or just represent a past, more energetic, period of activity. We 
also see part of an unusual cold front to the S. This region is seen 
clearly in the unsharp mask images (Fig. 2) and in one generated 
from data from all chips (Fig. 13). This structure appears to be con- 
nected to a region to the SW of the nucleus where the channel ap- 
pears. It could represent gas associated with subcluster merging in 
the cluster. Most likely given the relationship with the bubble chan- 
nel, the gas could be evidence of past energetic bubbles. The bubble 
channel is good evidence that the bubbles are not easily disrupted, 
prestmiably due to the magnetic structure (De Young et al 2003) 
and/or viscosity in the surrounding gas (Resmolds et al 2005). We 
assume that the pressure dips in the channel because there is unseen 
buoyant relativistic plasma there from the radio outbursts. 

An overall picture of the region is shown in the image of 
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Fig. 14, where data from all chips has been used. The structure 
of the inner regions can be seen together with the outer S bay, em- 
bedded within the more extended peak of cluster X-ray emission. 
The upper part of the Ha fountain (Fig. 2) can also be more clearly 
seen. 



7 SUMMARY 

Using a very deep, 900 ks Chandra image of the core of the Perseus 
cluster, we have found new outer features 50-80 kpc from the nu- 
cleus and measured the detailed properties of gas near the centre. 
The features are in the form of a concave cold front and baylike re- 
gion of hot gas which is in approximate pressure equilibrium. This 
could be the result of an energetic past outburst from the nucleus, 
or where bubbles accumulate. 

The inner radio bubbles are surrounded by complete higher 
pressure bands of gas behind a sharp front. The gas temperature 
does not change across the shock front, probably indicating that 
thermal conduction operates efficiently there, or that co-existing 
relativistic plasma mediates the shock. Pressure variations coinci- 
dent with ripples previously found in unsharp mask images reveal 
the presence of isothermal sound waves. The isothermal nature of 
the innermost gas means that a simple temperature estimate there 
does not reveal the expansion velocity of the bubbles. We suspect 
that they expand in rapid steps associated with outbursts of activity 
from the central radio source. Provided that the energy in the rip- 
ples is dissipated by viscosity, then the present heating rate in the 
ripples is sufficient to balance radiative cooling. Larger pressure 
variations are seen along a N-S channel, suggesting a sequence of 
bubbles, revealing the activity of the central radio source for the 
past 10^ yr. 

The gas in the centre is significantly multiphase with a large 
mass of gas 10^ Mq) associated with the optical Ha filamentary 
nebula, with ten times more mass in 0.5 keV gas than that radiating 
the optical emission lines. Mixing is likely occurring between the 
hot ICM and the cold filamentary gas, with much radiative cooling 
probably taking place below 10^ K. 

Cluster cores are complicated with the behaviour dependent 
on the bubbling of a central radio source and on microphysical 
transport processes. These in turn depend on the magnetic field 
structure, which itself may be a consequence of past cooling and 
bubbling. 
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Figure 14. Total 0.5-7 keV image from all Chandra CCD chips. 
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